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The absence of a missing moment inertia in clean solid 4 He suggests that the minimal experi- 
mentally relevant model is one in which disorder induces superfluidity in a bosonic lattice. To this 
end, we explore the relevance of the disordered Bose-Hubbard model in this context. We posit that 
a clean array 4 He atoms is a self-generated Mott insulator, that is, the 4 He atoms constitute the 
lattice as well as the 'charge carriers'. With this assumption, we are able to interpret the text- 
book defect-driven supersolids as excitations of either the lower or upper Hubbard bands. In the 
experiments at hand, disorder induces a closing of the Mott gap through the generation of mid-gap 
localized states at the chemical potential. Depending on the magnitude of the disorder, we find that 
the destruction of the Mott state takes place for d + z > 4 either through a Bose glass phase (strong 
disorder) or through a direct transition to a superfluid (weak disorder). For d + z < 4, disorder is 
always relevant. The critical value of the disorder that separates these two regimes is shown to be 
a function of the boson filling, interaction and the momentum cut off. We apply our work to the 
experimentally observed enhancement 3 He impurities has on the onset temperature for the missing 
moment of inertia. We find quantitative agreement with experimental trends. 



I. INTRODUCTION 

While superflow in a state of matter possessing a shear 
modulus might initially seem untenable, experimental 
claims for precisely this phenomenon in solid 4 He now 
abounci^^^^^ 1 ^ 1 ^ 3 .. Reported in the experi- 
ments by Kim and ChaniiS, (KC) was a dramatic change 
below 200mK in the period of a torsional oscillator con- 
taining solid 4 He. Because superfiuids come out of equi- 
librium and detach from the walls of the rotated con- 
tainer, they are expected to give rise to a period shift 
in such a geometry, assuming, of course, the rotation 
velocity is less than the critical velocity to create a vor- 
tex. The result is a missing moment of inerti a 14 ' 15 (MMI) 
and hence the period of oscillation decreases. The mag- 
nitude of the MMI is a direct measure of the superfluid 
fraction. In the original experiments reported by Kim 
and Chan 1 ^, the superfluid fraction ranged from .14% for 
4 He in vycor— to 2% in bulk 4 He. However, Rittner and 
Reppy 8 .^ have shown that the quench time for solidifying 
the liquid is pivotal in determining the superfluid frac- 
tion. In one extreme, when the sample is fully annealed, 
no MMI occurs. In the other, the MMI increased to 
an astounding 20% in samples in which the solidification 
from the liquid occurred in less than 2 minutes. While 
not all groups^ 3 , have been able to eliminate the MMI 
signal entirely by annealing 5 -!^ the sample and in fact 
there is at least one claim of MMI in a single crystal, 
the enhancement of MMI by a rapid quench does not 
seem to be in question. In fact, two independent experi- 
ments point to the key role played by disorder: 1) the To- 
doschenko et. al<i£ measurement that the melting curve 
of 4 He remained unchanged from the T 4 law expected 
for phonons in ultra-pure samples with a 3 He concen- 
tration of 0.3ppb and 2) the experiments of Clark and 
Chan 3 that increasing the 3 He impurity concentration 3 
from 20ppm to 40ppm increases the transition tempera- 
ture from 0.35K to 0.55K. 

Clearly the standard textbook supersolid in which 



vacancy or interstitial defects Bose condens ed 8 ' fails 
to explain the disorder dependence of the MMI. In 
fact, it is unclear at this writing even if a super- 
component is needed— to explain the MMI, primarily 
because experiments— designed to detect persistent mass 
flow have revealed no telltale signature. Monte Carlo 
simulations^ 2 - reveal, however, that superflow in solid 4 He 
is confined to grain boundaries. This observation is sup- 
ported by the experiments of Sasaki et alJ^ who observed 
mass flow only in samples containing grain boundaries. 
Nonetheless, the precise relationship between this exper- 
iment and the torsional oscillator measurements is un- 
clear because mass flow was observed at temperatures 
(1.1K which is not far from the bulk superfluid tran- 
sition temperature) vastly exceeding the onset temper- 
ature for MMI in the torsional oscillator experiments 2 , 
namely T c = 0.2K. 

Even if the MMI is not tied to superflow, disor- 
der is still the key player underlying the experimental 
observations^ 3 ". As disorder can induce superfluidity in 
the disordered Bose-Hubbard model, we explore its util- 
ity as a minimal model for the experimental observations. 
Certainly, this model does not have all of the microscopic 
details necessary to describe 4 He, in particular the pre- 
cise details needed to describe a grain boundary or the 
long-range interactions between 4 He atoms. Our central 
claim is only that it serves as a minimal model to de- 
scribe disorder-induced superflow in a bosonic system. 
Our work is based on a simple claim: 4 He is a hexag- 
onally close-packed self-generated Mott insulator. In a 
self-generated or self- assembled Mott insulator, the lat- 
tice and the 'charge carriers' are one and the same. In 
contrast, in Fermionic Mott insulators, the electrons oc- 
cupy pre-existing lattice sites formed by the ions. Our 
characterization of 4 He as a self-generated Mott insulator 
is relevant for three reasons: 1) In a supersolid the rel- 
evant transport is of the 4 He atoms themselves. Hence, 
if they form a Mott insulator in the clean system, no 
transport is possible. 2) Experiments^ 3 - and simulations 
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find an absence of MMI in the clean limi t 22 ' 23 . 3) We 
can immediately classify the candidate supersolids with 
this scheme because disorder can either— 1) self-dope 
the syste m 25 ' 26 or 2) create mid-gap states 2 ^. The former 
would generate either vacancies or interstitials and hence 
excitations in either the lower or upper Hubbard bands. 
The Andreev/Lifshitz 1 ^ scenario in which vacancies or 
interstitials Bose condense can be thought of as arising 
from doping a self-generated Mott insulator. We call such 
a state SSI. In electronic systems, disorder is well-known 
to have such an effect^ 5 .. We will show that SSI does not 
obtain in the disordered Bose Hubbard model. Rather a 
superfluid state (SS2) forms from overlapping localized 
mid-gap statea 27 i 28 i 29 > 30 i 31 i 32 i 33 . 34 i 35 . We argue that SS2 
is most relevant to the experimental observations. 

We establish several new results in this paper. First, 
we use the replica technique coupled with a renormalisa- 
tion group analysis to show that weak and large disor- 
der disrupt the Mott insulator (MI) in radically different 
ways. In particular, the critical value of the disorder 
that separates these two regimes is a decreasing func- 
tion of filling. Second, in the weak disorder regime a 
direct transition from the superfluid (SF) to the Mott 
insulator is possible whereas such a transition always in- 
volves the Bose glass (BG) phase at large disorder. This 
result resolves the controvers y 31 ! 33 ' 34 surrounding when 
the destruction of the superfluid necessitates an interme- 
diate Bose glass phase. Finally, we offer a quantitative 
test of this model by applying it to the 3 He enhance- 
ment of T c . The quantitative agreement suggests that 
the essence of the MMI in the experiments is captured 
by the disordered-Bose Hubbard model. 

II. INITIAL CONSIDERATIONS 

To describe boson motion in a random potential, we 
adopt the site-disordered Bose-Hubbard model 

H = -t^2 (b\bj + c.c) + + \ n i( n i - !)■ 

(1) 

In this model, b\ is the creation operator for a boson 
at site i and rii is the particle number operator and t 
and V are the Josephson coupling and on-site repulsion, 
respectively. We also define a value J — zt where z is the 
number of nearest neighboring sites. 

Though much of the theoretical work^iS&^i 3 ^! 3 ^ 5 . on 
the disordered Bose-Hubbard model has confirmed the 
originally proposed picture that an intermediate Bose 
glass localized phase disrupts the MI-SF transition, sev- 
eral key issues remain. 

PI. 7s there a direct MI-SF transition in the presence 
of disorder? 

For example, several analytical treatment s 27 ' 32 ' 34 ' 35 
suggest that the Bose-glass phase completely sur- 
rounds the Mott insulating phase, making a direct 



transition from the MI to SF impossible. How- 
ever, simulation s 28 ' 29 ' 31 and a renormalization group 
analysis 33 find that a Bose glass is absent in d = 2 at 
commensurate fillings. In fact, the renormalization group 
analysis of Pazmandi and Zimanyi 33 lays plain that the 
weak and strong disorder cases are fundamentally differ- 
ent. Only in the strong disordered case does the Bose 
glass phase completely surround the Mott lobes. How- 
ever, Herbut 3 ^ has also provided a convincing treatment 
of the large-filling limit and concluded that disorder is 
always relevant and destruction of the superfluid obtains 
through the Bose glass even in d = 2. 

P2. Do Mott insulators vanish for unbounded distri- 
butions? 

Fisher, et ali^l argued that no Mott insulating phases 
are possible when the width of the disorder exceeded V/2 
at T = 0. Consequently, for unbounded distributions, 
Mott insulators are absent at T = and only a super- 
fluid phase exists 27 -. Does the same hold for finite tem- 
perature? As the distributions characterizing disorder- 
in optical lattices are typically unbounded, this question 
must be resolved. 



A. Resolution 

We resolve both of the problems in this paper. 

First, we show that the missing ingredients that 
squares these seemingly contradictory results in PI are 
1) dimensionality, 2) critical momentum cutoff A c and 3) 
a filling and interaction-dependent critical value of the 
disorder A c . For e = 4 — (d + z) > 0, disorder is always 
relevant. In this case, the Mott insulating phase is de- 
stroyed and a BG obtains. This is in agreement with the 
work of Herbut^l on the destruction of superfluidity in 
d = 1 and 2 always takes place through the Bose glass. 
He finds that z = 1.93, implying that e > 0, matching 
our criterion for the relevance of disorder. 

For systems with e < 0, there exists a boundary in 
phase space separating disorder relevant and disorder ir- 
relevant regions. For filling m = 1, a direct transition is 
always allowed. For large fillings, the situation is more 
complicated. If the momentum cutoff (determined by 
the lattice constant) exceeds a critical value, A > A c , the 
MI is surrounded by a BG phase and direct transition 
from MI to SF is forbidden as illustrated in Fig. (TT]). In 
the opposite regime, A < A c , the strength of the disor- 
der is the key ingredient. For the weak disorder case, 
A < A c , a direct transition is allowed for large fillings 
while it is forbidden for strong disorder, A > A c . These 
results are in accord with the RG analysis of Pazmandi 
and Zimanyi 33 who studied an infinite range model and 
found^ 3 . that for e < a direct transition is possible. For 
e > 0, they found that disorder is in general relevant 
except at perhaps the particle-hole symmetric point at 
small filling where a direct transition survives at weak 
disorder. 

In addition, we analyse a Gaussian distribution for the 
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FIG. 1: Phase diagram for the disordered Bose- Hubbard 
model as a function of chemical potential fi/V and hopping 
strength, J/V. MI,BG and SF stands for Mott insulator, Bose 
glass and superfluid respectively. In the present of disorder, 
the lobes are shrunk, and we have two phase insider the lobes, 
MI and BG and outside, we have SF phase, a) The typical 
phase diagram when e < and A > A c or when e < 0, A < A c 
and A > A c . In this case, direct transition from MI to SF is 
only possible at m = 1. b) The typical phase diagram when 
e < 0,A < A c and A < A c . In this case, direct transition 
from MI to SF is possible for many filling numbers. 



site energies here and demonstrate how temperature and 
disorder are intertwined. At finite temperature, we estab- 
lish the existence of integer-filling Mott states. However, 
the T = analysis is beyond the scope of the treatment 
here as it corresponds to the infinite disorder limit. In 
particular, our replica analysis on unbounded distribu- 
tions is valid strictly when 



PA 2 /V < 1 



(2) 



where A is the variance of the distribution and V the 
on-site repulsion. In fact, this breakdown is fundamen- 
tally related to our central point that for bosonic systems, 
disorder destroys Mott insulators and gives rise to super- 
fluids. To see how this comes about, it is sufficient to 
integrate out the randomness by using the replica trick, 
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lim 

n->0 



Z n - 1 



(3) 



in which n represents the number of replicas and Z is 
the partition function. Performing the integral over the 
disorder 
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results in an effective Hamiltonian for the disordered 
problem. Here a is the replica index and we have assumed 
that the disorder is described by a Gaussian distribution 
of width A. We see clearly that the on-site interaction is 
replaced by 



V cS = V-f3A 2 



(•5) 



Consequently, at sufficiently low temperature, disorder 
can destroy the Mott gap. 

This paper is organized as follows. In the next sec- 
tion we compute the phase diagram for the disordered 
Bose Hubbard model using replicas and a renormaliza- 
tion group analysis. Since we start our analysis from the 
strongly interacting regime, any diagrams that are cal- 
culated cannot be computed using Wick's theorem. To 
circumvent this problem, we resorted to the analysis de- 
tailed in the Appendix. We explicitly compare the results 
for Gaussian and the uniform distribution case studied 
earlier 27 . The analysis of the phase boundary for the 
Bose glass is presented at the end of this section. This 
analysis is particularly lengthy as the topology of the 
phase boundaries is found to be delicately determined by 
the strength of the disorder and the cutoff. We close with 
an application of our central result that disorder enhances 
superfluidity to the problem of 3 He- induced enhancement 
of the onset temperature for missing moment of inertia. 



III. PHASE DIAGRAM OF THE DISORDERED 
BOSE HUBBARD MODEL 

In this section, we derive the phase diagram for the dis- 
ordered Bose-Hubbard model for the Gaussian and uni- 
form distribution of site energies. To establish the phase 
boundaries for the Mott insulator (MI) and superfluid 
(SF) phases, we employ a saddle-point analysis on the 
partition function^^ ^ 2 ., 

Z = Z J Y[D^(T)D^*(T)eM~S(^)} (6) 

i 

— In ( T T exp 



-nl>i(T)b: 



at 



R.c. 



(7) 



by introducing a Hubbard-Stratonovich field ipj to re- 
lease the b\bj term. Appearing in Eq. © are [J -1 ]^, 
the inverse matrix of hopping rates which will deter- 
mine the band structure for the kinetic energy and Zq = 



-Hot 



Trexp(-/?ffo), bf{r) = e H " T b^{0)e 

Differentiating the free energy with respect to tp yields 
the saddle point equation^ 



i 



(8) 



Because ipf is linearly related to (bf), its average value 
will serve to define the superfluid order parameter. This 
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can be seen more clearly by performing the cumulant 
expansion on &" (t). The action can then be rewritten 
as, 



+ « £ iciW+wi 4 

r ij = [J- 1 ] ij - 5aJ* j^drdr 1 (Tbl\r)b1{T')) (9) 

where r matrix acts as the mass term and hence deter- 
mines the appearance of superfiuid phase. 



A. Gaussian disorder 

For the Gaussian case, the Hamiltonian consist of two 
parts, 



tfo = ^f£«) 2 

a,i 
{i,3) 



£< 



(10) 



where /x e ff = M + V/2. Because the hopping term is a 
perturbation, our theory is valid strictly for V > J. In 
addition, since we are working in the limit in which the 
Mott lobes are well-formed, we must have that V e s > 
and PV 3> 1. The latter two constraints can be written 
as 1 > a where a — f3A 2 /V. It is this parameter that 
we will use to characterize the strength of the disorder. 
Using the eigenstates of H , that is, the eigenstates of 



particle number, (m\8) = ^e^^ m 6 ', we have 



Tbf(r)bW)) n = 







(11) 



x Y Z [{ m \e HaT K ] ~ e ~ Ho1 ~e HaT 'b b 3 e- H ° T '|m) 0{t - r') 

+ (m\e HaT 'b)e- HaT ' e^bfe-^^m^e^' - r)] . (12) 

For the above to be nonzero, we have to choose a = b 
and i = j. Inserting a complete set of states, 1 = 
ricSm<:=i \fn c )(m c \, between b^ and bf, we have only 
two terms left, \m a ± l(c = a)m c {c ^ a))(m a ± l(c = 
a)m c (c 7^ o)|. Note that we have replica symmetry 
between the initial and final states. However, replica 
symmetry breaking must be present in the intermediate 
states to have a nonzero correlation. The inserted state 
together with the creation and annihilation operators will 
lead to a term of the form Eo(mf ± 1, rrq) — Eo(mf, m\) 
where Eo(mf, rrq) is the eigenenergy of Hq. The explicit 



form for this term is 



(13) 



After integrating over r and t', we obtain, 

JdrJdr'(Tbf(r)b b J (r')) o 
(m + 1) ( 1 \ m / 1 



(m + 1) 



+ 



m 



(14) 



where we can neglect 1 /(3e± only when the temperature 
is small relative to the Mott gap, that is, ksT/e± -C 1. 
In the above equation, we are considering the energy of 
one replica, so the m°mj term will give rise to (n— l)mf , 
part of which is linear in n. Because we will take the 
limit n — > in the end, we can neglect all the high order 
terms when we calculate the energy of one replica. In 
terms of D = /3A 2 /2, the energies e± are defined as 



e± = E (m«±l,m*>)-E (m«,m*>) 

V eS (V cS r J& 2 \. 
= — ±ml — -(n-1)— J //<„'! 



(15) 
(16) 



e+(m) = mV — ji — (m + 1)1? 
e_(m) = (1 — m)V + fj,+ (m— 1)D (Gaussian) 



We defined m to be the integer closest to ii e s/V e s because 
in the low temperature limit, only this term in Z dom- 
inates. This holds for a system with non-conserved or 
commensurate particle number. For a system with con- 
served and incommensurate particle number, we should 
replace m by the particle number m, on each site. 

For the single-component case, r is a scalar and we just 
need r < to have superfiuid order. In our case, however, 
r is a matrix which must be diagonalized. For simplicity, 
we consider only nearest neighbor hopping in one dimen- 
sion case where Jy = t(5ij+i+5ij-i). The diagonal hop- 
ping matrix will be Sij Jcos( 2 ^) with j = 0, 1, N — 1 
which is the quantum number of momentum k = 2-jrn/L. 
Here N is the number of sites, L the system size, and 
J = zt where z is the number of nearest neighbors. Di- 
agonalizing the hopping matrix will of course require var- 
ious linear combinations of the ipi fields. Such linear com- 
binations will leave (rb^ (r)bf (t')^ invariant because of 

the appearing in front. Consequently, the condition 
for superfiuid order is, 



nj(n) = 



Note superfiuid order arises anytime one of the r\-s < 0. 
The phase diagrams we construct in this section corre- 
spond strictly to phases in which ipi = and Vi 7^ 0. 
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In the Bose glass section, we will make the distinction 
between the localized phase being gapped or ungapped. 

The onset of a MI state is determined by the largest 
eigenvalue of [J -1 ]. For a continuous band, this corre- 
sponds to 1/ J. Consequently, the phase boundary sepa- 
rating MI and SF phases is given by 



1 

J 



(to + 1) 



m 

E- 



(17) 



Using Eq. (TT5) . we rewrite Eq. (fTT)) as 

m(m - l)V 2 + V [(1 - 2m) n + J + D(m + 1 



2m 2 )} 



fj,[J + 2mD] - (to+ l)DJ+ (to 2 - l)D 2 = 



This equation describes a set of super-planes in terms of 
V — A — [i for different to. For a given chemical potential, 
it describes the phase boundary as a function of disorder 
and V. For to = 1, that is, one boson per site, we re- 
cover exactly 14(A) (Eq. ([36]) ) as the phase line between 
SF and MI. The analogous expressions can also be de- 
rived for fixed disorder a = D/V but varying chemical 
potential y = fx/V and x — J/V which reads, 



y 



X 

2 



ma 



± - - 2a) 2 + (4to + 2) (2a - l)x + x 2 



The result is shown in Fig. {2}. From the figures, we see 
that the distance between the upper and lower bound- 
aries of each lobe have shrunk by 2a and the whole lobe 
is shifted downward by ma relative to the ordered so- 
lution. As is evident, the MI phase still exists at finite 
temperature for the unbounded distribution. Finally, in- 
creasing disorder decreases the size of the Mott lobes. 
That the size of the Mott lobes shrinks with disorder has 
also been found in the extensive simulations of Trivedi 
and colleagues 2 ^ for a uniform distribution of site ener- 
gies. 



B. Uniform Distribution 

For completeness, we also compute the uniform distri- 
bution of site energies of width 2A studied in the original 
treatment of the disordered Bose-Hubbard problem 27 . 
Integrating over the disorder in this case is also straight- 
forward and yields 



n/ a - 

n / cm*!.-**- 

i 

-t J2 { b ? b i + - Z> + V/2)nt(19) 

(i,j},a i,a 

|E«) 2 -^ lnsinh (^ A E<) 
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FIG. 2: Phase diagram for the disordered Bose-Hubbard 
model with a Gaussian distribution of site energies. The two 
values of the disorder correspond to a = (3 A 2 /V = 0.4 and 
a = 0.6. 



H ,3 ?" : ') 



where the last two terms are interactions generated by 
(18) the integration over the disorder. Note the (3A 2 reduc- 
tion of the on-site repulsion is absent in the uniform dis- 
tribution case. Consequently, the T = limit can be 
taken explicitly. Introducing and still choosing the 
basis that diagonalizes Hq to perform the cumulant ex- 
pansion, we compute the last two terms at T — and 
the n — > limit to be 



lim — 



[i In sinh(/3A £ m? ± 1) - ~ ln(/3 £ m? ± 1)] 



[i In sinh(/3A £ m?) - ~ ln(/3 ^ m?)] 



— lim lim — In 
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sinh(/3AnTO^ ± /3A) 
sinh(/3Anm°) 
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lim — lncosh(±/3A) 

1 P 



P- 
= -A 

Thus we have, 



e+(m) = toU — /i — A 

e_(m) = (1 - m)V + n - A (uniform) (20) 

where \i is replaced by fi + A in the first term and by 
fi — A in the second term. It is this structure that makes 
the width of the MI lobes shrink by A as a function of 
filling relative to that in the clean limit. We then use Eq. 
(TH)) to obtain 
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FIG. 3: a) Phase diagram for disordered Bose-Hubbard model 
with uniform distribution with S — A/V = 0.4. b) Compar- 
ison between the replica theory and the treatment of Fisher, 
et al.— for the infinite-range hopping model with a uniform 
distribution of site energies. 



± -V(I - 25) 2 + (4m + 2) (28 - l)x + x 2 (21) 

as the phase boundary in the x — y(J/V — (J-/V) (x and y 
represent J/V and /i/V he) plane for the Mott insulator- 
superfluid transition. Here 8 = A/V. 

The phase diagram in the x — y plane shown in Fig. 
([3^,) bares close resemblance to the finite temperature 
counterpart of the Gaussian distribution. The only dif- 
ference between the two is that the disorder in the Gaus- 
sian case is characterized by a = /3A 2 , whereas for the 
uniform distribution at T = 0, the strength of the dis- 
order is set by 8 — A/V. Consequently, in the uniform 
distribution, the Mott lobes display a vertical shift of 8 
rather than a as in the Gaussian case. For an indepen- 
dent check on the accuracy of the replica method, we 
consider the uniform distribution but with infinite range 
hopping. In Fig. (J3t>) we compare the replica method 
with the mean-field criterion 
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>( m -l-(y+6)y 



(22) 



derived by Fisher, et al.— . As is evident, only minor 
quantitative differences obtain, lending credence to the 
replica treatment presented here. 



C. Bose Glass 

In the dirty boson model, a localized phase (Bose glass) 
exists in which disorder rather than the on-site repul- 
sion (Mott insulator) is the root cause. Unlike tradi- 
tional spin glass phases which are characterized by an 
Edwards- Anderson order parameter, the Bose glass does 



not admit such a description. In fact for the Bose- 
Hubbard model, the only Edwards- Anderson parame- 
ter that could be non-zero is (b1(t)b1(t')) . For the su- 
pcrfluid phase, this order parameter is trivially non- 
zero. However, there is no phase in which such order 
exists without simultaneously relying on superfluid or- 
der. With nearest-neighbour Coulomb interactions, such 
a glass is possible 43 , independent of superfluidity. The 
current analysis is limited, however, solely to the on-site 
Coulomb case. 

To analyze the Bose glass, we use the 
standar d 27 i 44 ' 45 i 46 one-loop renormalization group 
equations in conjunction with the mean-field phase 
boundaries to derive a criterion for the onset of the Bose 
glass phase. The field theory of our model is, 
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(23) 



where ao is the lattice constant. The coefficients g aa and 
g a b can be calculated using the cumulant expansion pro- 
cedure outlined in the Appendix. For a Gaussian distri- 
bution, these coefficients are given by 
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(24) 
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4 [(m + |)F - (m + 2)D - n] 
m(m — 1) 



(m- §)V + (to-2)D + /x 



(25) 



The signature of the disorder-induced localized phase is 
the divergence of the coupling constant for the interaction 
between different replicas. To this end, we derive the one- 
loop renormalization equations 46 



dg a 



dg a b 



(flu 



gac9ca I 



eg a b + K d ((4 + 2p)(g aa + g b b)g a b 



gaegea I 



for the coupling constants g a b and g a b- Here £ is the 
standard rescaling parameter, e = 4 — (d + z), K e i = 



7 



(47r)d /2 r(d/2) , (K 2 = P is the number of the compo- 
nent of ip ( p = 2 in this case) and d is the spatial dimen- 
sion. We are particularly interested in p = 2 and d = 2. 
Care must be taken in analyzing these equations, how- 
ever, as the coupling constants, g a b and g aa are actually 
ultrametric matrices. Using the Parish multiplication 
rule for such matrices, we partition g a b into a diagonal 
part g and an off-diagonal part which is a function g(x) 
defined in the domain x G (0, 1). At the replica symmet- 
ric fixed point, we find that 
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ep 



9{x) = - 



16(p-l)K d 
e(4-p) 



16(p - l)K d 



(26) 

for i£ [0,1]. (27) 



This fixed point is unstable" forp > 4(1 — e). That is, for 
d < 3.5, there is a runaway to the strong disorder region 
signalled by g(x) — > oo, the signature of localization. The 
main criterion for the boundary to separate the disorder 
relevant and disorder irrelevant regions comes from the 
renormalization equation for g{x). If we consider the 
replica symmetric case, we only need two parameters, 
the off-diagonal, g(x) = g, and diagonal parts, g. The 
renormalization equations simplify to 



dg 
dg_ 



= eg-K d [(p + 2)g 2 +pg 2 ] 

= {e-K d (4 + 2p)g}g + (4-2p)g 2 . (28) 



Two characteristic properties of the Bose glass in this 
RG scheme are 1) g — > oo, and 2) ip = 0. As is well 
known, when g — ► +oo, the RG procedure breaks down. 
Hence, we can use the RG procedure to demarcate the 
boundary between the disorder relevant and disorder ir- 
relevant regimes. From Eq. (|28p . the condition for g to 
run to infinity is 



e - K d (A + 2p)g > 0. 



(29) 



So the boundary separating disorder relevant region and 
disorder irrelevant region is, 



e - K d (A + 2p)g = 0. 



(30) 



This result is in fact similar to the long wavelength 
limit derived by Fisher et al^I. In fact, they applied 
the replica trick and RG analysis to a similar mean field 
action. Without considering the p— dependence, they 
found that the coefficient of g 2 (see Eq. (|2"B)> ) is always 
positive. However, as is clear from Eq. (|28p. in the gen- 
eral case when the p— dependence is considered, this coef- 
ficient can be negative. Note the presence of p is two-fold 
as it also generates a cross term gg in the RG equations. 

Eq. (3D]) together with that for ip = 0, that is, r > 
give rise to the Bose glass phase boundary in the phase di- 
agram. This criterion depends on e, x — J/V, y = [i/V , 
disorder and the momentum cutoff. Hence, if e, the dis- 
order, and the momentum cutoff are given, Eq. (|30p will 



define a series of curves in the x — y plane. For Gaussian 
disorder, in the domain n/V G (to + i — 28, m — | + 8), 
g < 0. As a result, if e > 0, Eq. (f29|) is always satisfied, 
which means that for systems with d + z < 4, disorder 
is always relevant. In this case, the rp = regions all 
turn into the Bose glass and a direct transition between 
a Mott insulator and the superfluid is not possible. If 
e < 0, the criterion (Eq. (J23)) will separate disorder rele- 
vant and irrelevant regions in the x — y plane. In general, 
the criterion depends on x, y, the disorder, 5 and the 
momentum cutoff A, and is given by 



x d = x d (y,A,e) = 



-48n 2 e 



A 2 /V 

(m + l)(m + 2) 



I [m — (m + l)a — y] [m + | — (m + 2)a — y] 
m(m — 1) 

[1 - to + (to - 1)q + y] [§ - to + (to - 2)a + y] J 

For different fillings to, we have a class of curves which 
form concentric lobes x — x d (y, A, 5) in the x — y plane. 
So for each filling number m with e < 0, we have a critical 
value x d . If x > x d , disorder relevant, while for x < x d 
disorder is irrelevant. 

The tips of the MI lobes are at y = (to + l)a — l + (l — 
2a)\Jm(m + 1), precisely where x reaches its maximal 
value. Recall a = (3A 2 /V. For large fillings, m — > oo 
and y approach to yo where yo — m — \ — ma. To 
see whether the disorder-relevant region lies within the 
MI lobes in x — y plane (recalling that x = J/V and 
y = (J./V) , we calculate the value of x, x Ml , for the MI 
lobes evaluated at yo and x BG of BG lobes evaluated at 

BG 

yo. Consequently, we consider the ratio, ttt- This ratio 



-1/2 



„BG 



.MI 



2 7 r A /6(-e)(l - 2a) 



2m + 1 



TO + 1 



(31) 



is a decreasing function of filling number to, and for large 
filling number, 



„BG 



4 7 r- v /6(-e)(l - 2a) 



as m 



(32) 



Whether we can have a direct transition from MI to SF 
depends on the whether the above ratio is greater or less 
than one. If the above ratio is greater than one, it means 
that the curve demarcating the disorder-relevant region 
intersects the r = curve. In this case, the BG region is 
located in the upper and lower regions of MI lobes as de- 
picted in Fig. |T]). Consequently, in such cases, a direct 
transition from MI to SF is allowed. If the above ratio 
(Eq.(|3"2"])) is less than one, the BG surrounds the MI and 
a direct transition between the MI and SF is forbidden. 
Because for m = 1, x c runs to infinity as y — > 0, the dis- 
order curve always intersects the r = curves at to = 1. 
Consequently, we reach the conclusion that for m = 1, 
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a direct transition is always allowed. This prediction is 
in principle testable by direct numerical simulation. The 
phase diagram, Fig. in the A/J—V/J plane displays 
the direct transition from the MI to the superfluid as the 
disorder is increased. In this plane, a further increase in 
the disorder leads to a transition to m — 2 Mott insulat- 
ing state. Hence, we predict that the superfluid density 
should be a non-monotonic function of the disorder. Sim- 
ilar conclusions have been reached in a Landau-Ginzburg 
treatment^ of the supersolid problem. 

The ratio in Eq. (|3ip depends on the momentum cutoff 
A and the disorder a. We can see that increasing disorder 
a will decrease the ratio, so if A is less than a critical 
value, no matter what strength the disorder is, the ratio is 
always less than one; thus direct transitions are forbidden 
except for m = 1. So for a given momentum cut-off A 
and interaction V , a direction transition from MI to SF 
is forbidden if A < A c where 



A 



96n 2 (-e)V 2 



(33) 



which follows from < 1 assuming a = 0. In this case, 
for any disorder strength, a direct transition is impossible 
between the MI and the SF at large fillings. This corrob- 
orates the result derived earlier by Herbut^i that the MI 
phase is always surrounded by the BG in the large-filling 
limit of the Bosc-Hubbard model. 

BG 

For large momentum cutoff A > A c , the ratio ^mt 
could be greater or less than one depending on the disor- 
der strength a. A critical value of A c exists. Hence, for 
weak disorder, A < A c , we have a direct transition for 
large fillings A < A c with 



£Aj 
V 



1 



f 



yl 



967r 2 (-e) 

which is derived from < 1 assuming A > A c . 
Renormalization also modifies the value of r t< 



(34) 



,46 



r(0 = r exp{[2-K d ((2+p)g+pg(x)]Z} 

= r cxp j-^J =r i?- 1/2 . (35) 

where we have only considered the replica symmetric 
case. Here R c is the correlation length R c = exp(£) oc 
r^ v w r Q 1 , so r(£) oc r^ 4 which means that the 
renormalization does not shift the MI-SF phase boundary 
which occurs at r = 0. 

Ultimately, it is the Bose glass that makes the disor- 
dered boson problem distinct from the disordered elec- 
tronic Mott insulator. In the presence of disorder, the 
boson lattice adjusts (contracts or expands) so that the 
chemical potential remains in the gap. In the electron 
problem, in which the electrons occupy pre-existing lat- 
tice sites, disorder changes the position of the chemical 
potential 25 -. Consequently, for the boson problem, it is 
the nature of the in-gap states that ultimately determines 
whether the disordered system is localized or not. How- 
ever, as we see here the criterion is a complicated function 
of the system parameters. 



4« 



Ml (m=2) 




Ml (m=1) 



10 
V/J 



15 



20 



FIG. 4: A typical phase diagram in A/ J — V/ J plane for 
disordered Bose-Hubbard with Gaussian distributed disorder. 
There are three phase: MI,BG and SF. Phase at "X" which 
is a MI becomes a SF when the disorder is increased. 



D. He impurities 

3 Hc increase the onset temperature for the missing mo- 
ment of inertia. Although we do not have a microscopic 
model for a grain boundary, the point defect model we 
have outlayed here explains this effect qualitatively as 
disorder can enhance the superfluid region. In essence, a 
disordered system with interaction V can be represented 
by a pure system with an effective interaction V c g. If 
for a pure system, the critical interaction is V c , then for 
a disordered system, the corresponding critical point is 
Vcs — V — j3A 2 = V c . An immediate consequence is that 
the new boundary for the Mott-superfluid transition is 
shifted to higher values of the on-site interaction. That 
is, for the disordered system (with one boson per site) V c 
is replaced by 



V e (A) = V C + p c A 2 



(36) 



Consequently, to remain on the phase boundary, increas- 
ing the disorder must be compensated by an increase 
in the onset temperature as is seen experimentally 3 - for 
3 He defects and studied theoretically by a Abrahams and 
Balatsky^ using a Landau-Ginzburg approach. To for- 
malise this, we consider 3 He defects with a concentration 
c and an on-site energy e^. We will treat the 4 He atoms as 
having on-site energy e± with concentration 1— c. A rigor- 
ous treatment require a binomial distribution of disorder. 
But to get the basic scene of the influence of disorder, we 
use Gaussian distribution to approach this disorder. The 
key parameter is the variance of the distribution of on- 
site energies, A 2 = A 2 , + c(l — c)(e 2 — £i) 2 , where A 2 , is 
the disorder which can be eliminated by annealing. For 
a clean system, the transition from the Mott insulator to 
the superfluid is given by k B T c /J = (V c - V)/V<M-. We 
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FIG. 5: Critical temperature as a function of impurity con- 
centration. Experimental data are taken from Ref. (7). 



now replace V by V e s and solve for T c . The solution, 



K B T C = piJ 



+ P2J + P3 Jc(l- c), (37) 



has a square-root dependence with pi = 



and p 3 



_ CVc-V) 



2V C 



= yf 



(g2-gl) 



-y . Knowing that the critical temperature 

is 0.2K in the absence of 3 He impurities leads to a rela- 
tionship between P2J and p\J . Thus, we have two free 
parameters p\ J and p 3 J to fit the curve. We show in Fig. 
(|5|) a plot with the fitting parameters: p\ J = — 0.10K, 
p 3 J = (90K) 2 and p 2 J = (0.28K) 2 . From the above for- 
mula, we can see that if there are no impurities and no 
other disorder that can be annealed away, that is, both 
c = and = 0, we obtain a negative T c which means 
there is no supersolid transition. Also, if disorder is too 
large, V e s = V — /3A 2 < 0, and the "net interaction" is 
attractive which results in an insulating phase. Conse- 
quently, for sufficiently large disorder, we also obtain an 
absence of a supersolid transition. Hence, although the 
treatment here is not rigorous, it sufficiently rich to cap- 
ture the interplay between disorder, finite temperature, 
and supersolidity. The quantitative agreement, which is 
tied more to the functional form than the fitting param- 
eters, lends credence to our claim that disorder underlies 
the missing moment of inertia in solid 4 He. 



IV. CONCLUSION 

We have presented what we think is the minimal model 
that captures disorder-induced superfluidity in bosonic 
systems. While we undoubtedly do not have the suf- 
ficient microscopic details to model actual grain bound- 
aries, the results presented here offer a general framework 
in which the general problem of disorder-induced super- 
fluidity can be formulated consistently. We have seen 
from our replica analysis and the one-loop renormaliza- 
tion analysis that the phase boundaries of the disordered 
Bose-Ffubbard model can be determined but do not ap- 
pear to be universal, in contrast to the phase boundaries 



constructed from general considerations in the early work 
of Fisher, et al^i. In particular, a direct MI-SF transition 
is possible as found earlie r 31 i 33 ; however, the criterion 
depends on the disorder, interaction strength and filling 
numbers. Further, we have shown how Mott insulating 
phases can be observed in unbounded distributions. This 
application is particularly relevant to experiments 3 ^ on 
optical lattices as the disorder in such systems always 
obeys an unbounded distribution. Since the flows are 
to the strong disorder limit, a treatment (currently not 
available) in this parameter space is essential to under- 
standing the phase structure of the Bose-Hubbard model. 
Finally, because the MI phases always give rise to super- 
fluids for intermediate disorder (for example, < D < V 
for Gaussian distributions), we believe this model is the 
correct starting point for analysing the reports of miss- 
ing moment of inertia in solid He induced by disorder, 
in particular the extreme sensitivity of the critical tem- 
perature to 3 He impurities. 
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V. APPENDIX 

In the section, we will derive the effective action. The 
action is, 



i.a^b 



(38) 



where in momentum space, 
1 



Jcos(fcao) 



- / dT(T T b(T)b\0)) (39) 



and ao is the lattice constant. The coefficients u,v are 
given by the averages, 

U = ~J dT 1 ...dT 4 (T T b a (T 1 )b a (T 2 )b a \T 3 )b a \T 4 )) 



24 



dT 1 d7a(r T 6 a (ri)6°t(7j I )) 
rfT 1 ...dr 4 <T T 6 Q (T 1 )6 b (r 2 )6 a t( r3 )6ft (T4) ) 

dndr 2 (T T b a ( n )b^(T 2 )) 

dT 3 dT 4 (T T b b (T 3 )b b ^(T 4 )) 



(40) 



(41) 
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To compute these correlation functions, we insert a com- 
plete set of states 



IliJm? ><mf| = 1 



(42) 
(43) 



in between all b ± (ri) operators and integrate over all Tj. 
The terms from the first term of the order of u, v of which 
the order of replica indices are aabb or bbaa will cancel 
with the second term. So we have, 



24 



(m + l)(m + 2) m(m — 1) 
£i(ei + e 2 ) 



(m+ l) 2 



m(m + 1) 



6e !,o( e i.o + £ 2,i) 12e 2 (e^ + e ,i) 
m(m + 1) 

12e_i )0 (e_i,o + eo,-i)(e_i,o + eo,-i - eo,i) 

to(to + 1) 
12ei i o(ei,o + eo,i)(ei,o + £ o,i — e o,-i) 
to(to + 1) m 2 



12e 2 _ 1)0 (e_i,o + eo,-i) 6e^ X)0 (e_i,o + £-2,-i) 
where the energies are defined as follows: 

E{m1,m\) 



^(mf + TO f)- Meff ( m ? +m f) 



/3A 2 



(44) 



eo,-i 
ei,o 

£2,1 

e-2,-i 

e -l,o 

eo,i 

e±i 
e±2 



BK + l^J-lJ-EK.m^-l) 
- (m + 2)D - 

mY — (to + 1)_D — /i 

£(m? + 1, mj + 1) - £(r<, mf + 1) 

mV — mD — /i 

E{m1 - 1, to^ - 1) - E(mt, m\ - 1) 

— (to — 1)V + mD + /i 

£(m? - 1, mj) - E{m%, mj) = e_ 

-(to - 1)Y + (to - 1)L> 

E{m a % - 1, + 1) - E(m a t) m\ + 1) 

-(m- 1)V + (m- 2)D + n 

E{m1 ±l,TO 4 b ) - £(m?,m?) = e± 

£(m? ± 2, mj) - #(ro?±, m\) 



e_i + e_2 



ei + £2 



-("»- -)y+(m-2)£) + /i 
(m+i)V- (m + 2)D-n 



(45) 



Here we have used the trick lim„^o Sf>^ a A = 
lim„^o(n— 1)^4 = —A We can simplify the result further 
to 

(to + 1) 2 to 2 



12e\(e++ D/2) I2et fs_ + D/2) 
m(m + 1) ( 1 I ' 



12(F-3£>) ' £-7 (46) 

To proceed, we make the approximation 

l/[Jcos(fca<>)] « 1/[J(1 - (fca ) 2 /2)] = (1/J)[1 + {ka f /2\A7) 

We then rescale ^ by the factor ip — > y/Ja^ ^tp an d 
replace the J^-, /3 by / / d/3 respectively. We have, 



Sty) = J ' d d x J ' d/3 { x; 



1 - J 



m + 1 to 



2 , 



x |^ a (^r)| 2 + i|V^(x,r)| 2 + ^ E |^ a (x,r)| 



J 2 v 



Xl^ a (x,r)| 2 |^(x,r)| 2 + 0(|^| 6 )l (48) 



a^6 



Denoting the momentum cutoff A = — . from above we 
have, 



J 2 



J 2 A 2 



9ab = —oU 



I2n 2 

m(m + 1) / 1 



(m + l) 2 



m 



e\{e++D/2) e 2 _(e-+D/2) 



(V-3D) \e 



S- 



(49) 



,/ 2 



J 2 A 2 
"48tt 2 



(m + l)(m + 2) 



4 [(m + |)F - (to + 2)£> - n] 



m(m — 1) 



e 2 . [-(to - §)V + (to - 2)L> + ^] 
where the dependence on the cut-off is allowed. 
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